Measuring the Quantum State 
of a Large Angular Momentum 
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We demonstrate a general method to measure the quantum state of an angular momentum of 
arbitrary magnitude. The (2F + 1) x (2F + 1) density matrix is completely determined from a set 
of Stern-Gerlach measurements with (4F + 1) different orientations of the quantization axis. We 
implement the protocol for laser cooled Cesium atoms in the 6Si/2(F — 4) hyperfine ground state 
and apply it to a variety of test states prepared by optical pumping and Larmor precession. A 
comparison of input and measured states shows typical reconstruction fidelities T > 0.95. 
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The quantum state of a physical system encodes in- 
formation which can be used to predict the outcome of 
measurements. The inverse problem was mentioned by 
Pauli already in 1933 Q : is it possible to uniquely deter- 
mine an unknown quantum state by measuring a suffi- 
ciently complete set of observables on a number of iden- 
tically prepared copies of the system? This very basic 
question has gained new relevance in recent years, fol- 
lowing the realization that systems whose components 
and evolution are manifestly quantum can perform tasks 
that are impossible with classical devices, such as certain 
computations, secure communication and teleportation 
1^ . As we harness quantum coherent dynamics for such 
purposes, the development of techniques to accurately 
control and measure quantum states becomes a matter 
of practical as well as fundamental interest. Reconstruc- 
tion of a (generally mixed) quantum state based on a 
record of measurements is a nontrivial problem with no 
general solution ||, but system-specific algorithms have 
been developed and demonstrated experimentally in a 
limited number of cases. These include light fields [@, 
molecular vibrations |)j , electron orbital motion Q , and 
center-of-mass motion in ion traps |?J and atomic beams 
[^j . More recently multi-particle quantum states have 
been measured for entangled spin- 1/2 systems in NMR 
H and for polarization-entangled photon pairs [|l0| . 

In this letter we present a new experimental method 
to measure the unknown quantum state for an angu- 
lar momentum of arbitrary magnitude. The protocol 
is implemented for laser cooled Cesium atoms in the 
6£i/2 (F — 4) hyperfine ground state, and typically re- 
produces input test states with a fidelity better than 0.95. 
Our work is motivated in part by our ongoing study of 
quantum transport and quantum coherence in magneto- 
optical lattices, where the atomic spin degrees of freedom 
couple to the center-of-mass motion [[y]]. The correla- 
tion between spin and motion in this system offers the 
possibility to use the angular momentum as a 'meter' in 
the sense introduced by von Neumann [ ft2"[ , to probe the 
spinor wavepacket dynamics. A similar application has 



been proposed in cavity QED, where an atomic angular 
momentum can be used to read out the quantum state of 
light in an optical cavity pH| . We expect our technique 
also to provide a powerful experimental tool to evaluate 
the performance and error modes of quantum logic gates 
for neutral atoms ifuj . 

The angular momentum quantum state of an ensemble 
of atoms with spin quantum number F is described by its 
density matrix g. Newton and Young have shown that 
sufficient information to determine g can be extracted 
from a set of 4F + 1 Stern-Gerlach measurements, car- 
ried out with different orientations of the quantization 
axis |]l5|| . They derived an explicit solution for quantiza- 
tion directions with a fixed polar angle 6 and evenly 
distributed azimuthal angles Our reconstruction al- 
gorithm uses the same general approach, but takes ad- 
vantage of a simple numerical method to solve for the 
density matrix for a less restrictive choice of directions. 
This allows for flexibility in the experimental setup and 
improves the robustness against errors. As a starting 
point for our reconstruction we choose a space fixed co- 
ordinate system {x, y, z} in which to determine g. If we 
perform a Stern-Gerlach measurement with the quantiza- 
tion axis along z, we obtain the populations of the 2F+ 1 
eigenstates \m z ) of F z , i.e., the diagonal elements of the 
density matrix g. Information about the off-diagonal ele- 
ments can be extracted from additional Stern-Gerlach 
measurements with quantization axes along directions 

fik z. For each of these measurements we obtain a new 

(k) 

set of populations ir m , corresponding to the diagonal el- 
ements of a matrix g^ representing g in a rotated coor- 
dinate system with the quantization axis oriented along 
fik- The associated coordinate transformation is a rota- 
tion by an angle 9k around an axis Uk in the xy-plane and 
perpendicular to hk [see Fig. 0(a)]. The new populations 
can then be found from a unitary transformation of g, 
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where —F < m,i,j < F and 1 < k < 4F + 1. The 
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FIG. 1. (a) Direction hk of a Stern-Gerlach measurement 
in spherical coordinates, (b) Time-of-fiight distribution with 
well separated peaks corresponding to the populations in the 
magnetic sublevels \m). 



rotation operators are given by = exp — iOkFiik 

Arranging all the populations 7T™ and the density ma- 
trix elements into vectors n and p, Eq. [l] can compactly 
be written as ir = Mp, where the elements of the rect- 
angular matrix M are determined by the set of rotations 
The elements of the density matrix g can then 
be obtained as p = M + 7r, where 



(2) 



is the Moore-Penrose pseudoinverse [jl6| . R and Xi are the 
rank and non-zero eigenvalues of the Hermitian matrix 
M^M, with Wi and v$ being the corresponding eigenvec- 
tors of M^M and MM 1 , respectively. 

In a real experiment there is both random measure- 
ment noise and systematic errors, which cause the ob- 

(k) 

served populations ir) n ' to deviate from those predicted 
by Eq. [l| In this situation the pseudoinverse solution 
yields a least-squares fit to the data |L6]]. This fit, how- 
ever, will become sensitive to noise and errors if one or 
more of the Xi in Eq. || are too close to zero. A robust 
reconstruction algorithm must therefore use a set of di- 
rections {hk} that lead to reasonably large A^. A second 
problem arises because the pseudoinverse solution, while 
always the best fit to the data, is not guaranteed to be a 
physically valid density matrix when noise and errors are 
present — a density matrix has unit trace, is Hermitian 
and has non-negative eigenvalues. The pseudoinverse so- 
lution automatically fulfills the first two conditions when 
the input populations are normalized. However, when 
negative eigenvalues occur we have to employ a different 
method to solve the inversion problem. For this purpose, 
we decompose g as 



q = TT f with Tij = 0, for j > i 
and Tu £ 



Ei r «i a = 1 - 



This form automatically enforces the density matrix to 
have unit trace, be Hermitian and positive semi-definite. 
The (2F + l) 2 — 1 independent real- valued parameters 
of the complex, lower-triangular matrix T can then be 



optimized to yield the best fit between measured and 
expected populations. 

We have implemented this general procedure for Ce- 
sium atoms in the 6Si/2{F = 4) hyperfine ground state 
manifold. A standard magneto-optical trap (MOT) and 
optical molasses setup is used to prepare an ensemble 
of ^10 6 atoms in a volume of 0.1 mm 3 and at a tem- 
perature of 3.5/dK. Atoms from the molasses are loaded 
into a near-resonance optical lattice, composed of a pair 
of laser beams with orthogonal linear polarizations (ID 
lin_Llin configuration) and counter-propagating along the 
(vertical) z-axis. Following this second laser cooling step 
the atoms are released from the lattice, at which point a 
range of angular momentum quantum states can be pro- 
duced as discussed below. To allow accurate preparation 
and manipulation of these test states, we measure and 
compensate the background magnetic field in our setup 
to better than 1 mG. 

A laser cooling setup provides a convenient framework 
for Stern-Gerlach measurements |l7|]. At the beginning 
of each measurement we define the quantization axis by 
applying (switching time ~2 ps) a homogeneous bias 
magnetic field of ~1 G pointing in the desired direc- 
tion. All subsequent changes in the magnetic field are 
adiabatic, which ensures that the initial projection of 
the atomic spin onto the direction of the field is pre- 
served at later times. As the atoms fall under the in- 
fluence of gravity, we apply a strong magnetic field gra- 
dient (|B| « 100 G and V|B| « 100 G/cm ) by pulsing 
on the MOT coils for 15 ms. The resulting state depen- 
dent force F = —mgppB^7\^\ is sufficient to completely 
separate the arrival times for atoms in different \m) as 
they fall through a probe beam located 6.9 cm below 
the MOT volume [see Fig. [l](b)]. The magnetic popu- 
lations can then be accurately determined from a fit to 
each separate arrival distribution. In our case a total of 
4F +1 = 17 Stern-Gerlach measurements are needed to 
reconstruct g. The corresponding 17 different directions 
hk of the bias magnetic field are produced by three or- 
thogonal pairs of coils in ncar-Hclmholtz configuration. 
This setup provides complete freedom to choose the mea- 
surement directions, and allows us to set the dc- magnetic 
field with an accuracy of 0.1° in a volume of 1 cm 3 . For 
the results reported below we use 16 measurements on 
a cone with 9^ w 82°, and a single measurement with 
0k = 0°- 

We evaluate the performance of our reconstruction pro- 
cedure by applying it to a number of known input states. 
These test states are created by a combination of laser 
cooling, optical pumping and Larmor precession in an 
externally applied magnetic field. First consider opti- 
cal pumping by a er~-polarized laser beam tuned to the 
651/2 (F = 4) -> 6P 3/2 (F' = 4) transition, which al- 
lows us to prepare the ensemble in a pure \m z = —4) 
state. A single Stern-Gerlach measurement with quan- 
tization axis along z confirms the essentially unit pop- 
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FIG. 2. (a) Density matrix of a pure \m z = —4) input state 
obtained by optical pumping with (j _ -polarized light, and (b) 
measured density matrix with a fidelity of T = 0.97. Note: 
All figures display absolute values of the density matrices. 



ulation of this state. For a single non-zero population 
we can rule out off-diagonal elements due to the con- 
straint \gij\ 2 < QuQjj, and in this case we therefore know 
the complete density matrix with a high degree of con- 
fidence. Both the input and measured density matrices 
are shown in Fig. ||, and show excellent agreement. To 
quantify the performance we use the fidelity 18 



T 



Tr 



1/2 1/2 
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which is a measure of the closeness between the input 
(f?i) and the reconstructed (g r ) density matrices and 
takes on the value T — 1 when they are identical. For 
an \m z = —4) input state our reconstruction fidelity is 
T = 0.97. Optical pumping with a + yields a nearly 
pure \m z — 4) state, which we can similarly reconstruct 
with a fidelity of T = 0.94. Starting from \m z = —4) 
we can further produce a range of spin-coherent states 
p9[ by applying a magnetic field along, e.g., the x-axis 
and letting the state precess for a fraction of a Larmor 
period. Figure |^ shows the measured density matrices 
for four different precession angles, again with excellent 
reconstruction fidelity. 

It is desirable to also check the reconstruction of test 
states that are not quasi-classical, and whose density ma- 
trix exhibits large coherences far from the diagonal. One 
example of such a state is \m v = 0), which we can pro- 
duce by optical pumping on the F = 4 — > F' = 4 tran- 
sition with linear polarization along y. In this basis, the 
input test state has the density matrix shown in Fig. |](a). 
Its representation in the {Iro^)} basis is easily found by a 
(numerical) rotation by 90° around x [Fig. |](c)], an d can 
be directly compared to the reconstruction [Fig. ^(d)]. 
To simplify the visual comparison we finally rotate the 
measured density matrix by —90° around x to find its 
representation in the {Im^)} basis [Fig. |](b)]. The ba- 
sis independent reconstruction fidelity is T — 0.96. The 
non-classical nature of the measured state is apparent in 
its Wigner function representation pQ |, which takes on 
negative values as shown in Fig. |^. 

As a final test we have applied our procedure to the 




FIG. 3. Measured density matrices for spin-coherent states 
obtained by a controlled Larmor precession of |m = —4) 
around x by (a) 30° \T = 0.95], (b) 60° \T = 0.96], (c) 90° 
\T = 0.95], and (d) 120° \T = 0.92]. 



(presumably) mixed states that result from laser cool- 
ing. Figures §|(a,b) show the input and reconstructed 
states produced by laser cooling in a ID lin_Llin lattice. 
The input state shown here is based on a single Stern- 
Gerlach measurement of the magnetic populations, but 
it seems reasonable to assume that the highly dissipative 
laser cooling process destroys coherences between mag- 
netic sublevels and that the density matrix is diagonal. 
Our measurement confirms this assumption. We can per- 
form the same experiment with atoms released directly 
from a 3D optical molasses, in which case there is no pre- 
ferred spatial direction and one expects something close 
to a maximally mixed state. Figures ^(c,d) show input 
and measured density matrices in good agreement with 
this assumption. 

In summary we have demonstrated a method to ex- 
perimentally determine the complete density matrix of 
a large angular momentum, and tested its performance 
with a variety of angular momentum quantum states 
of an ensemble of laser cooled Cesium atoms in the 
6S1/2 (F = 4) hyperfine ground state. The input and 
reconstructed density matrices typically agree with a fi- 
delity of J- > 0.95. Limitations on the accuracy to which 
the density matrix can be measured appear to derive from 
ss 3% uncertainties in the individual population measure- 
ments, and from variations in the direction of the bias 
magnetic field during the first few ^s when it is turned 
on. These variations are likely due to induced currents 
in metallic fixtures and additional coils in the vicinity 
of our glass vacuum cell, and translate into uncertainty 
about the exact orientation of the quantization axes for 
the Stern-Gerlach measurements. Efforts are underway 
to eliminate these problems and achieve even better re- 
construction fidelities. 
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FIG. 4. (a) Input state \m y = 0) prepared by optical pump- 
ing with linear j/-polarized light, and the reconstruction result 
(b) represented as density matrices in the {|m w )} basis. The 
density matrices of input (c) and measured state (d) in the 
standard {|m z }} basis. The fidelity is 0.96. 
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FIG. 5. (a) Measured \m y = 0) state, represented by the 
Wigner function W(0, <p) in spherical phase-space. The darker 
shading indicates negative values. Selected cuts of W(8, tp) 
along (b) (p = —90° and (c) 6 — 90° (dotted lines = input 
state, solid lines = measured state). 
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